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Interacting quantum systems with strong or long-range interactions 
exhibit quite a rich phase structure. Cluster formation and supercon- 
ductivity are examples. These phenomena are also expected in the astro- 
physical context, e.g., during the formation or in the structure of neutron 
stars. To describe these phenomena a proper treatment has to go be- 
yond the simple picture of noninteracting quasiparticles. An appealing 
formalism for a systematic approach is provided by the framework of 
Zmi. "Dyson equations. Within an equal (imaginary) time formalism Dyson 
^ equations can be derived for an arbitrary large cluster embedded in a 
medium Q|. For practical use and the sake of simplicity the medium 
£N) is treated as uncorrelated to derive the respective n-body cluster Green 
q functions. Further, we neglect "backward" propagating particles, so the 
qj Fock spaces for different number of particles n are disconnected. This 
way it is possible to derive effective in-medium n-body equations that 
tan be solved rigorously with few-body techniques 0|l3| ■ 
. These resulting two-, three-, and four-body equations elaborated here 
include the dominant medium effects in a systematic way. These are 
I the self energy corrections for masses and the Pauli blocking that in 
J> turn leads to a change of binding energies, viz. change of the masses of 
f — flusters, and change of reaction rates. Further, within this approach the 
OO critical temperatures for condensation (of bosons containing two or four 
particles) are calculated. 

Defining Ho — YH=i £ i with the quasi-particlc self energy 
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G (z) = (z - Ho)- 1 N = R (z) N. (2) 



ei = fc?/2mi + 2V2(12,12)/2 (1) 
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and the Fermi function fi = /(ei) = l/(e /3 ( £l ~ Al ) + 1), the n-particle 
cluster resolvent Go is 
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Here Go, Ho, and N are matrices in n particle space and z denotes the 
Matsubara frequency |14[ . The Pauli-blocking factors for n-particles are 



N=fJ 2 ...f n ±f 1 f 2 ...f n , f = l- f (3) 

Note: NR = R N. Defining the effective potential V = £ pairs a N 2 V 2 
the full and the channel resolvents are 

G{z) = (z-H - = R(z)N, (4) 

G a (z) = {z-H - N^)- X N = R a (z)N, (5) 

Note that V ] j= V and R(z)N ^ NR{z). For the scattering problem it 
is convenient to define the in-medium AGS operator Up a (z) 

R{z) = S a pR (z) + Rp{z)U Pa (z)R a (z) (6) 

that after some algebra leads to the in-medium AGS equation 

U Pa (z) = SpMz)- 1 + J2 S fJl N^(z)R (z)U ia (z), (7) 
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where 5p a = 1 — 5p a . The square of this ^-operator is directly linked 
to the differential cross section for the scattering process a — ► (3. The 
driving kernel consists of the two-body t-matrix derived in the same for- 
malism, however given earlier and known as Feynman-Galitskii equation 



A numerical solution using a coupled Yamaguchi potential has been given 
in Ref. For a temperature T = 10 MeV and the three-body system at 
rest in the medium results a given in Fig.^ For the bound state problem 
it is convenient to introduce form factors 

l^)=E^WIfe>. (8) 
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Since the potential is nonsymmetric right and left eigenvectors are dif- 
ferent, although the bound state energies are the same, 

\F a ) = hfiNilf(B 3 )Ro(B 3 )\Fp), (9) 

\F a ) = 5 a 0Tg(B 3 )NPRo(B 3 )\Fp). (10) 
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FIG. 1. Break-up cross section for different densities of nuclear matter for 
temperature T = 10 MeV 



The binding energy depends on fi, T, P c . m .. Results for P c . m , = 
are given in Fig.^J for different potentials and temperatures. Note that 
the dependence on density is rather similar for two different potentials 
studied, although the binding energies for the isolated triton differs by 
10% the Mott density is practically at the same place once the binding 
energies are renormalized to each other. For helium the Mott density 
is smaller due to the Coulomb force, however for asymmetric nuclear 
matter, e.g. NJN n ~ 0.72 (for the 129 Xe+ 119 Sn reaction) this effect is 
compensated 115]. The dependence of the Mott effect on the momentum 
is given in Ref. [ 5| . 

We now turn to the four-body problem in matter. In addition to hav- 
ing different channels as for the three body system now the channels 
appear in different partitions that makes the four-body problem even 
more involved. The partitions of the four-body clusters are denoted 
by p, r,a,.. ., e.g., p = (123) (4), (234) (1), ... for 3 + 1-type partitions, 
or p = (12)(43), (23) (41), ... for 2 + 2-type partitions. The two-body 
sub-channels are denoted by pair indices a,/3, 7,..., e.g. pairs (12), 
(24),. . . The two- and three-body i-matrices have to be defined with re- 
spect to the partitions that leads to additional indices. The four-body 
in- medium homogeneous AGS equation are defined for the form factors 
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FIG. 2. Triton binding energy as a function of nuclear matter density. 
Dashed-dot-dot T = 10 MeV, other T = 20 MeV. 



=E^EW" 1 (B4)|^>, (11) 

r a 

where 5p a — Sp a , if (3, a C r and (5^ Q = otherwise and \ip a ) is the 
a-particle wave function. They read ^ 

\ T D = J2^tU^(b 4 )r (b 4 )n2t2(b 4 )r q (b 4 )\t;), (12) 

where a C er, 7 C r. A numerical solution of this equation is rather 
complex. In order to reduce computational time we introduce a en- 
ergy dependent pole expansion (EDPE) that has been proven useful in 
many application involving the a-particle and is accurate enough for the 
present purpose. However, we have to generalize the original version of 
the EDPE because of different right and left eigenvectors. Details will 
be omitted here, see || 

In the two-body sub-system the EDPE reads 

T l( Z ) - ^2 \r~fn(z))t~fn(z)(r 7n (z)\ ~ ^ \9~tn)hn(z){g ln \ 
n n 

= ^N2\g ln )t rn {z)(g ln \. (13) 
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where we have chosen a Yamaguchi ansatz for the form factors for sim- 
plicity. Inserting this ansatz into the Feynman-Galitskii equation deter- 
mines the propagator t jn (z). In the three-body sub-system the EDPE 
expansion reads 

(g 0m (z)\Mz)U^(z)Mz)\9,n(z)) * £ \rZ(z)Kl(z)(T^(z)\. 

(14) 

with the three-body EDPE functions 

|f$*(*)> = (^^0(^)1^)^(53)1^), (15) 

that we get from solving the proper Sturmian equations 

^|f«/> = ^2{9a n \Ro(B 3 )\gp m )tp m {B 3 )\T^) (16) 

/3m 

m^r^T) = ^{g a MB 3 )\gp m )t Pm {B 3 )\r^) (i?) 

Inserting everything into the homogeneous AGS equations allows us to 
redefine the form factors that are now operators in the coordinates of 
the 2 + 2 or 3 + 1 system, respectively 

K S ) = E(r^ m , I ,(S4)|^ m (S4)(5^n(S4)|i 1 , 0(S4)|^} (18) 

0m 

and therefore the final homogeneous equation 

n s ) = EEE^ r (r-r(54)ii 7 „(S4)if^(54)) o^) i r ?>> 

rt uk •yn 

(19) 

is an effective one-body equation with and effective potential V and an 
effective resolvent Q - 

vp rt W = J2^r(r^(z)\t in (z)\r^(z)), (20) 

g;i%\z)=W- (2i) 
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FIG. 3. Energy dependence of the binding energy of the a-particle 

The binding energies of the two-, three-, and four-body systems are 
shown in Figjsj for a temperature of T = 10 MeV and a cm. momentum 
of the respective cluster of P c . m . = 0. The B = line reflects the respec- 
tive continuum. Investigating the zeros of the two-body Joost function 
the quasi deuteron survives as an anti-bound state (not resonance) with 
increasing densities, viz. for energies above the continuum flq| . The fate 
of the triton and of the a particle for B > still needs to investigated 
as well as a possible appearance of Efimov states related to B — > of 
the sub-system. Since the Efimov states are 'excited' states, e.g. for the 
three-body system close to the 2 + 1 threshold, their blocking may be 
smaller since the wave functions contain higher momentum components. 
Note that the slope of the binding energies as a function of densities 
for the larger clusters is also larger. This is a clear indication that the 
masses of the clusters change with increasing density. The Mott density 
of the a particle appears at 

B±{n Mott = 0.19fm- 3 ,T = 10 MeV,P c . m . = 0) = 

For comparison a perturbative result using Gaussian functions fitted to 
the charge radius of the a-particle is also given. 
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FIG. 4. Temperature vs. effective chemical potential (/i* = (x — S(0)), 
E(0) Hartree-Fock shift of the single particle energy at zero momentum. Lines 
show the critical temperature for pairing (solid) and quartetting. The dashed 
line shows a result given in JTt| , the diamonds show the solution of the AGS 
equation given in ( |19[ ) 

Finally, I address the question of a possible four-particle condensate 
or quartetting 0. The condition is S 4 (n, T c , P c . m . = 0) = 4/x. From 
FigJ^ we argue that a condensation is likely, i.e. 

where the critical temperature for a condensation turns out to be higher 
than for the pairing. However, for /i > the situation seems not so clear, 
since the four-body AGS equation ( jl9| ) develops poles related to zeros in 
4/i — H . Unlike the two-body case were these poles disappear because 
the numerator becomes as well zero at 2 fx and T c (viz. 1 — /j — / 2 —> 0). A 
vanishing numerator is not obvious for the four-body case because there 
are more channels involved. It remains to be clarified, if the rapid fall 
of the critical temperature found in Ref. |fl7| using variational treatment 
with square-integrable functions remains, if one uses an exact treatment 
of the four-body problem. This is currently investigated. 
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